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The three dimensional mean spherical model on a hypercubic lattice with a film geometry Lxoo 
under periodic boundary conditions is considered in the presence of an external magnetic field H . 
5-H , The universal Casimir amplitude A and the Binder's cumulant ratio B are calculated exactly and 

found to be A = -2C(3)/(5ir) « -0.153051 and B = 2tt/( % /5 ln 3 [(l + \/5)/2]). A discussion on 
the relations between the finite temperature C-function, usually defined for quantum systems, and 
the excess free energy (due to the finite-size contributions to the free energy of the system) scaling 
' function is presented. It is demonstrated that the C-function of the model equals 4/5 at the bulk 

critical temperature T c . It is analytically shown that the excess free energy is a monotonically 
increasing function of the temperature T and of the magnetic field \H\ in the vicinity of T c . This 
' property is supposed to hold for any classical d-dimensional 0(n), n > 2, model with a film geometry 

under periodic boundary conditions when d < 3. An analytical evidence is also presented to confirm 
that the Casimir force in the system is negative both below and in the vicinity of the bulk critical 
temperature T c . 
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I. INTRODUCTION 
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I The Casimir effect is a phenomenon common to all systems characterized by fluctuating quantities satisfying some 
conditions on the boundaries of the system (for a general review on the Casimir effect see, e.g., §, §). In the 
statistical mechanical systems the Casimir force is usually characterized by the excess free energy 

O ' 

/» b (T, L) = / a , b (T, L) - Lf hulk (T), (1) 

due to the finite size contributions to the free energy of finite systems with a film geometry L x oo 2 , where boundary 

• conditions a and b are imposed on the surfaces bounding the system across the direction L. Here fa. b(T, L) is the full 
lO ! free energy per unit area (and per ksT ) of such a system and /bulk is the corresponding bulk free energy density. 
,— I • The Casimir force 

m ■ 

OC, ./ Casimir U > lj ) ~ \ Z > 

• then arises naturally in the thermodynamics of these confined systems. 

For 0(n)-symmetric model systems (n > 1 ), depending on the boundary condition (a, b) and on n, f®\(T, L) may 
or may not contain contributions independent of L. For the Ising-like systems, i.e. n = 1, these can be the surface 
free energies f s ,a,(T) and f s ,b(T), and the interface free energy fi(T) (for brevity we consider the dependence on the 
temperature T only). For the 0(n),n > 2, models these will be only the contributions stemming from the surface 
free energies because the analog of the interface free energy is the helicity modulus T(T) and the corresponding 
contribution is of the order T(T)/L. In general at the critical temperature T c (of the corresponding bulk, i.e. L = oo, 
!LJ ' system) the full free energy f a ^(T,L) has the asymptotic form 



X ■ /a.b(T c , L) - Lf hulk (T c ) + fsATc) + f s .h(T c ) + A a . b + • • • , (3) 

where d is the dimensionality of the considered system and A aib is the so-called Casimir amplitude. The L-dependence 
of the Casimir term (the last one in Eq. (0)) follows from the scale invariance of the free energy and has been derived 
by Fisher and de Gennes B. The amplitude A a ,b is universal, depending on the bulk universality class and the 
universality classes of the boundary conditions 0], In the present article we will only consider the case of periodic 
boundary conditions (which implies that / Sia = / Si b = T = ). Then, according to the standard finite-size scaling 
theory (see, e.g., || for a general review), near the critical temperature T c and in the presence of a small external 
magnetic filed h the behavior of / ex is given by 

r(t, h, L) = L-^X^iattL 1 ^, a h hL A /»), (4) 
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where t — (T — T c )/T c is the reduced temperature, h = H/(ksT), at and ah are nonuniversal scaling factors, A cx 
is universal (usually geometry dependent) scaling function, X ox (0,0) = A pcr , and v and A are the corresponding 
(universal) scaling exponents. 

An interesting point of view on the properties of the excess free energy comes from the finite-temperature general- 
izations of the Zamolodchikov's C-theorcm j(| for quantum systems with arbitrary dimensionality due to Netto and 
Fradkin (7) (see also Zabzin ||; for a general review on phase transitions in quantum system see, e.g., flQ|| ). They 
define from the free energy a function C of the coupling constants and the temperature which is a positive and, in 
the regimes where the quantum fluctuations dominate, a monotonically increasing function of the temperature. The 
C function is, in fact, an analog of the excess free energy of the system that they consider. 

Before passing to a discussion of some details it seems necessary to comment on the well known point that for 
temperature driven phase transitions with T c > the quantum fluctuations are unimportant near the temperature 
critical point. Therefore, it seems that the properties of the system around one quantum critical point (with respect 
to a given quantum parameter, say, g) at T = cannot tell us anything about the properties of this system around 
its temperature critical point T c > 0. In fact the dimensional crossover rule helps to make a bridge between these 
phenomena. According to this rule the critical singularities (with respect to g, T — 0) of a rf-dimensional quantum 
system are formally equivalent to those of a d + z classical one (z is the dynamical critical exponent) and critical 
temperature T c > 0. On that idea are actually based the investigations of the low-temperature effects in quantum 
systems (see, e.g., P^tl), i.e. one considers an effective system with d infinite space and z finite "temperature" 
("imaginary-time") dimensions Lt ~ {Ti/ {ksT)) 1 ^ with periodic boundary conditions, and applies the methods of 
the finite-size scaling theory (in what follows we will set U = k,B = 1)- An exact lattice realization of these ideas is 
presented in fll5| . 

Since the generalizations of the Zamolodchikov's C-theorem are formulated for quantum systems with z = 1, in the 
remainder we will focus our attention on such class of systems only. For these systems Netto and Fradkin define Q| 
the dimensionless function 

C (ft g, a) = -fi d+1 n{d) lim V' 1 [F v (/?, g, a) - E (g, a)} , (5) 

where Eq is the zero-temperature energy, i.e. the energy of the "infinite" system, V is the volume (V — ► oo, but N/V 
is fixed, where N is the number of particles), n(d) is a positive real quantity, (3 — 1/T , Fy is the full free energy of 
the "finite" system (where the only "finite" dimension is the "temperature" one, i.e. the "geometry" of the system is 
oo d x Lt) and a is the characteristic length scale of the lattice. The real positive quantity h(d) is supposed to be of the 
form v d /n(d), where n(d) is a positive real number (which depends only on the dimensionality of the system) and v is 
the characteristic velocity (e.g. the velocity of the quasiparticles) in the system. Obviously, the exact choice of n(d) 
does not effect the monotonicity properties of the C-function. In JjJ the definitions n(d) = T((d+1) /2)((d+l) /n( d+1 ^ 2 
for bosons and n(d) — T((d+ l)/2)((d + 1)(2 — 2 1 ~ d )/n( d+1 ^ 2 for fermions have been proposed. 

In accordance with the dimensional crossover rule the statement that C is a positive and a monotonically increasing 
function of the temperature can be "translated" in a statement that the function — A cx of the corresponding classical 
system is positive and a monotonically increasing function of L _1 , see Eqs. (^) and (||) (of course, the last is equivalent 
to a statement, that A cx is a negative and a monotonically increasing function of L ). In R], Q it is shown that 
the monotonicity of the C-function is related to the absence of long range order in the systems under consideration. 
The existence of long range order destroys the general validity of the monotonicity. Within the classical systems no 
long-range order exists above their bulk critical point. So, we expect the statement formulated for A cx to be generally 
valid above T c for any classical system. Supposing that this is true and recalling that in the vicinity of T c X cx is a 
function of the scaling variables x\ = attL x l v and X2 = ahhL A ^ u , which both are monotonically increasing functions 
of i, we come to the conclusion that in the vicinity of its critical temperature (T > T c ) the excess free energy 
of a given system is a monotonically increasing function of any of its scaling parameters when the other one is kept 
fixed. Since x\ and Xi are monotonically increasing functions of the temperature and the magnetic field, respectively, 
the last implies that A cx , in the vicinity of T c , is a monotonically increasing function of t (t > 0) and h too. It 
is possible to present some arguments to support that the above statement can be extended to the region t < for 
0(n), n > 2, systems in contrast with the Ising-like systems. The reasoning for the difference in the expected behavior 
of the excess free energy in 0(n) and Ising-type models is closely related to the well know differences in the behavior 
of the correlation length £,oo{T) in these models: in the Ising model £,oo(T) < oo both below and above the bulk critical 
temperature, whereas in 0(n), n > 2, models below T c and in the absence of an external filed (h = 0), due to the 
existence of soft modes in the system (spin waves), £oc(T) is identically infinite. On that basis one expects that, away 
from T c , A cx will tend to zero exponentially fast in L (see, e.g., ||) for the Ising-type models, and, therefore, being 
of the order of around T c , X cx cannot be a monotonic function of its scaling parameters in the vicinity of T c . 

In 0(n), n > 2, models the finite size corrections should be essential not only in the vicinity but also below T c Q. In 
other words, we expect the monotonicity in the behavior of the correlation length in 0{n) 1 n > 2, models around T c 
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to be mirrored by a corresponding monotonic behavior of the excess free energy. If an external field is applied (h 7^ 0) 
then £oo(T, h) < 00 and, of course, we expect that X cx — > exponentially fast with L again, similarly to the Ising-like 
systems behavior. But, since X° x < 0, for any fixed t < the last implies that X CK will be a monotonically increasing 
function of the magnetic field in the under critical vicinity of T c too. 

The statements presented above, should be considered, of course, only as a plausible hypothesis, which has to be 
checked in order to probe the region of its validity. For example, it is under question if the monotonicity property of 
X CK will still hold if the finite system undergoes a phase transitions of its own. It is reasonable to believe that the 
hypothesis holds for any 0(n), n>2, system with d < 3 (then in the finite system with short range interaction there 
will be no "real" phase transition). 

In the present article we will show, within the three-dimensional mean spherical model, that in the vicinity of 
T c the excess free energy scaling function X°* is, indeed, a monotonically increasing function of any of its scaling 
parameters (x\ and x 2 ) when the other one is kept fixed. The last implies that X cyi is a monotonically increasing 
function of t, h and L above T c , and monotonically increasing, with respect to t and h, but monotonically decreasing, 
with respect to L, function below T c . 

Let us turn now to the behavior of /casimir- From Eqs. (||) and (Q) it immediately follows || 

/Casimir(i, h, L) = L~ d X C&simir (xi, X 2 ), (6) 

where the Casimir force scaling function is 

Id Ad 

^Casimir (x l,X 2 ) = (d - l)X^( Xl ,X 2 ) ~ ~X X — X eX (xi , X 2 ) X 2 —X^{x 1 ,X 2 ). (7) 

v ax\ v ox 2 

Note, that Xcasimir is again a universal function of x\ and x 2 . We remind that for finite-size systems this means that 
^Casimir will be the same for all systems of the same universality class and geometry and boundary conditions. It is 
believed that if a = b the Casimir force will be negative (see, e.g., jl7|l; strictl y sp eaking, for an Ising-like system 
this is supposed to be true above the wetting transition temperature T w [ft.6|, p7| , In the case of a fluid confined 

between identical walls this implies that then the net force between the plates will be attractive for large separations. 
One of the goals of the present article is to prove analytically this general expectation, i.e. that ^CasimirC^i, x 2 ) < 
for any [x\, x 2 ) € R 2 , on the example of one exactly solvable model. We will also show that if T < T c and H — the 
Casimir force is a monotonically increasing function of the temperature. We believe that these properties are valid 
for any 0(n), n > 2 models. 

The full temperature dependence of the Casimir force has been investigated exactly in two-dimensional Ising strips 
by Evans and Stecki whereas the upper critical temperature dependence of the force in 0(n) systems has 

been considered by Krech and Dietrich p(| by means of the field-theoretical renormalization group theory in 4 — e 
dimensions. (For the Ising-like case they have derived also some results for T < T c .) The only example where exact 
expression for the Casimir force as a function of both the temperature and the magnetic field is available is that one 
of the three-dimensional mean spherical model [^) . By numerical evaluation of the expressions derived there it has 
been shown that the force is negative, i.e. it is consistent with an attraction of the plates confining the system. The 
most results available at the moment are for the Casimir amplitudes A a ^. For two-dimensional systems at T — T c by 
using conformal-field theory methods the amplitudes are exactly known for a large class of two-dimensional models , 
|24| , (2^] . In addition to the "flat geometries" recently some results about the Casimir amplitudes between spherical 
particles in a critical fluid have been derived too |^6| . For d = 3 the results for the Casimir amplitudes available in the 
Ising-like case have been obtained by Migdal-Kadanoff renormalization-group calculations , by some interpolation 
of the exact values for d = 2 and d = 4 |^] , and, relatively recently, by Monte Carlo methods |^2| , [^3| . For n > 2 
the only existing results are obtained by the e-expansion technique, where the calculations are performed up to the 
first order in e pCfl . 

In the present article the hypotheses for the monotonicity of the excess free energy and that the Casimir force is 
negative under periodic boundary conditions will be verified analytically on the example of the three-dimensional 
mean spherical model. We will present also simple analytical results for the universal values of the Casimir amplitude 
and the Binder's cumulant ratio. If one takes the normalization factor of the analog of the C-function in the form for 
bosons (this will keep the C-function of the critical Gaussian model to be C — 1 for any d), it will be shown that the 
"C-function of the three dimensional spherical model" is 4/5 at the critical point. As it is well known, the infinite 
translational invariant spherical model is equivalent to the n — > 00 limit of the corresponding n-component system 

The results we are going to present are an extension and continuation of those published in H. In the notations 
and the definitions in the remainder we will closely follow ||. That is why here we only briefly recall, in Section 2, 
the definition of the model and give the final expressions, obtained there, for the excess free energy and the Casimir 
force which will be our starting expressions for the aims of the current article. In Section 3 we verify the hypotheses, 



3 



formulated above, for the excess free energy and the Casimir force. In Section 4 we derive the exact universal values 
for the Casimir amplitude and the Binder's cumulant ratio. The paper closes with concluding remarks given in Section 
5. 



II. THE MODEL 



We consider the ferromagnetic mean-spherical model (see, e.g., p8| , p9[ for a general review) on a fully finite 
e?-dimensional hypercubic lattice of | A | sites and with block geometry L\ x L 2 X ■ ■ • x Ld, where Li, i = 1, — , d 
are measured in units of the lattice spacing. The Hamiltonian has the form 



b.c. 



(WieA) = -\ K E jh it a i a i + s E°i~ 



(8) 



i£A 



= Jjf-), and h 



Here <Ji £ M, i € Ad (crj = a (r^)) is a variable, describing the spin on lattice site i (at r^), s is the spherical field, K 
is a dimensionless coupling, jf- c - is a matrix with dimensionless elements, so that (K/f3) Jfj c - is the exchange energy 
between the nearest neighbors (under boundary conditions b.c.) spins at sites i and j (of course, jfj 
is the external magnetic field. The dependence on the boundary conditions is denoted by a superscript b.c 

The scaling function of the free energy density of the spherical model has been discussed in details in the literature 
for different boundary conditions, dimensionalities and geometries of the system, for both the cases of short as well 
as for long range interactions in the Hamiltonian Q, p9fl , pfj| , pi] , |3^ ). By any of the approaches used there one 
can, of course, derive an expression for the excess free energy scaling function. Here, for d = 3 and under periodic 
boundary conditions we will take it in the form given in |j| 



1 



X^(x u x 2 ) = -(4tt) 



-3/2 



E 



(-l) fc frE +1 -l& +1 ) 

(k + 1)! (k - 1/2) 



(9) 



/OO 
dxx' 2 (1 + 2R(4tt 2 x)) cxp[~y L x] 



-2 / dxx 5/2 R(l/Ax)exp[-y L x] + j dxx 5/2 exp {—yoo%) 
Jo 

1 2 ( 1 1 \ 1 / 

"0^2 + o- Tl ~ ^ ■ 



where 



R ( x ) = E exp [~ x,?2 ] ' 

9=1 



(10) 



.1-1 



= (if- K c )L,x 2 = K- 1/2 hL 5/2 (11) 



are the scaling variables (note the difference in the definitions of X\ here and in the Introduction; now x\ decreases 
when T increases), 



I* 00 

K c = dx [exp (-2x) 7 (2a;)] 3 = 0.25273 . 
Jo 



(12) 



is the critical coupling, and yh and j/00 are the solutions of the spherical field equations that follow from (^j) by 
requiring the first partial derivatives of the right-hand side of ([)]) with respect to j/l and j/oo to be zero. 

For the finite-size scaling function of the Casimir force one immediately obtains from Eqs. (^), (^), (|^) and the 
definitions of the scaling variables x\ and x 2 J5| 

Acasimir (£1,3:2) = 2A cx (xi , x 2 ) - \x\ ( — — ] - \x x (y^ - y L ) . (13) 

2 VJ/oo VlJ 2 

The equations @ - (|l^) provide the basis of our further analysis. 
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III. VERIFICATION OF THE HYPOTHESES 



We will prove analytically that the finite-size scaling function of the excess free energy, given by Eq. (^|), is a 
monotonically increasing function of any of its scaling parameters X\ and x 2 when the other one is kept fixed. First, 
by using the identity 



E 

k=0 



(-i)V +1 

(ft + 1)! (k - 1/2) 



3 y 3 W 



,-3/2 



exp (— x) dx — -(1 — exp(— y)), 



the Jacobi identity for the R function (see Eq. (pL0| 



1 



R[A^x) =- 



> A~KX 



1 + 2R{ ~x 



- 1 



and taking into account that 

f°° dx / 1 \ 

J Q ^5/2 R ( ) ex p(-^) = A ^Wv Li 2 (exp (-Vv)) + Li 3 ( ex P {-Vv))] > 

after some elementary manipulations we obtain from (|^) 

1 



(14) 



(15) 



(16) 



X ex (x u x 2 ) = - 



2tt 

+\ x l [— ) +r;xi(y 00 -y L ), 



y L Li 2 (exp (-y/yZ)) + Li 3 (exp (-y/yT)) 



VL 



(17) 



where Li p (z) are the polylogarithm functions. The main advantage of the above representation of X ex is the existence 
of some nontrivial identities |33|], [Q for the polylogarithm functions (see next Section) that allow the universal 
constant A = X ex (0,0) to be expressed in a simple closed form. 

The spherical field equations for y^ and y^ can be now rewritten in the well known and very simple forms (see, 
e.g., for h = 0, Eq. (86) in M) 



xi_J_ 



2 sinh 



1 



and 



x\ 



47T 



(18) 



(19) 



where the first equation is for the finite and the second one for the infinite system, respectively. In order to obtain 
Eq. ( |l8| ) use has been made of the facts that dLi p (x)/dx = Li p _i(x)/a; and Lii(x) = — ln(l — x). Let us denote by 
9l(x2, yL) the right-hand side of Eq. ( ]l8] ) and by 500(^2,2/00) the right-hand side of Eq. (|l9|). Then, it is easy to see, 
that 



9L{x 2 ,y) = goc(x2,y) - - eacp(— y/y)] 



(20) 



From the above equation and having in mind that in Eqs. ( |18| ) and (19) j/l > 0, y > we conclude that 

9L{x 2 ,y) > goc(x2, y). (21) 

It is also elementary to verify that gL{x 2 ,yL) and floofe, S7oo) are monotonically decreasing functions of yL and ?/oo, 
respectively. Let now yoo{xi,X2) be the solution of Eq. ( |19| ) for given x\ and x 2 - Then, from E q. (|2l|) , the fact that 
giy(x2, y) is a monotonically decreasing function of y, and that for the solution yi,{xi, x 2 ) of Eq. (|18|) one should have 
gh{x 2 ,y L ) = goo(x 2 ,yoo), we obtain 

Vl(xi,x 2 ) > yoo(xi,x 2 ). (22) 

We are now ready to deal with the monotonicity properties of the excess free energy scaling function. From ( |l7| ) and 
having in mind the spherical filed equations ( [iq ) and (|l^) we derive 
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dX cx 
dxi 



-2 iVL 



and 



dX c 



1 

VL 



(23) 



(24) 



From these expressions and Eq. (|22|), taking into account the definitions of the scaling variables (11), we obtain that 
the excess free energy scaling function is a monotonically increasing function of both the temperature T and the 
magnetic field \H\. As a function of the finite size L of the system the scaling function is monotonically increasing 
above and decreasing below T c . These properties of the scaling function as a function of the scaling variables x\ and 
x 2 are illustrated in Fig 1. One clearly sees that for any fixed X2 the scaling function is a monotonically decreasing 
function of X\, and, for any fixed x\ a monotonically increasing function of (o^l- Finally, it is worth mentioning, that, 
for x 2 = from ( |l7| ) and ( p2| ) it immediately follows that X cx < 0. From Fig. 1 one observes that this is true also for 
x 2 + 0. 

We turn now to properties of the Casimir force. Our aim is to show that the force is negative under periodic 
boundary conditions for any values of T and H. The finite-size behavior of the Casimir force in the vicinity of the 
critical point is given by Eq. (^J) where the scaling function is given by Eq. @. For T < T c the same expressions 
are actually valid with the only difference that the definition of the variable x 2 now should be x 2 = K^^hL 5 ' 2 and 
x\ ^> 1. Here we are not going to discuss if then the above expressions can be simplified further e.g., being a function 
of a given combination of x\ and x 2 , as it is usually the case of first order phase transitions [pq| ). For T away above 
T c the Casimir force, as it has been shown in |5j, tends to zero exponentially fast with L in a full accordance with 
the general expectations about its behavior above the critical point. We will not be interested in the explicit form of 
these exponentially small corrections. Having in mind all these comments, for the behavior of the Casimir force for 
any T and H one obtains explicitly 



/Casimir(i, h, L) — L 



1 1 

VL Vac 



3/2 3/2 



Vo 



(25) 



+ \fyl Li 2 (exp (-- v/yT)) + Li 3 (exp {-y/yZ)) 



Since the inequality ( p2[ ) is still valid, from the above expression it immediately follows that fcasimh(t,h) < 0. 
Numerical evaluation of the behavior of the finite-size scaling function of the Casimir force has been given in j5j . It 
is in full agreement with our analytical result. Finally we show, that for T < T c and h — 0, i.e. Xi > and x 2 = 0, 
the Casimir force is a monotonically increasing function of the temperature, i.e. monotonically decreasing function of 
x\ . From (pq) and taking into account that j/oo = when T < T c we obtain 



dx\ 



a 



r (xi,0) 



1 



1 dy L 



dx\ 



(26) 



From ( ^8[ ) it is easy to see that dyh/dxi < 0, and, therefore c/Xcasiroir(^i, Q)/dxi < 0, i.e. the Casimir force is an 
increasing function of T for T < T c and h = 0. 

In this way we have completely verified the hypotheses formulated in the introductory part of the article for the 
behavior of the excess free energy and the Casimir force. 



IV. CASIMIR AMPLITUDE, C-FUNCTION AND BINDER'S CUMULANT RATIO 



Here we will be interested in the properties of the system at its bulk critical point. This implies xi = x 2 — with 
a solution of the spherical field equations (see Eqs. (|n]) and fll8|)) j/oo = and vl = yL, c = 4 In 2 [(l + /2] (this 
value of yL.c is well known and, seems, has been derived for the first time in |p6[). The problem of determination of 
the Casimir amplitude reduces now to exact evaluation of the expression 



•(CO) --JL 



(exp ( 



(exp (- 



(27) 



Denoting by r the "golden mean", i.e. r = (l + v5j /2, it is easy to show that 



G 



exp {-y/yZ^) = T 2 = 2 - r, 
which reduces the above problem for X cx (0,0) to the problem for evaluation of the expression 



(28) 



Li 3 (2 - r) - In (2 - r) Li 2 (2 - r) - - In 3 (2 - r) 



(29) 



Fortunately, this is exactly the problem solved by Sachdev 34 studying his example of a conformal field theory in 
three dimensions. By the help of some polylogarithm identities he has shown that a — 4£(3)/5. Therefore, we obtain 
for the Casimir amplitude of the three dimensional spherical model under periodic boundary conditions 



2C(3) 

5-7T 



-0.153051. 



(30) 



The numerical value of this amplitude has already been reported in ||. Recalling now that — JT cx (0, 0) corresponds 
to the analog of the C-function for our model and taking the normalization factor in the form that will keep the 
C-function of the critical Gaussian model to be C = 1 for any d (i.e. by taking the normalization in the form for 
bosons) we conclude that the "C-function of the spherical model" is 4/5 (at T = T c for d = 3 under periodic boundary 
conditions). 

Let us turn now to determination of the Binder's cumulant ratio for the considered model. We will use for it the 
definition of the form fa] (up to a prefactor 1/3) 



Br = -L 



_ d X {4) (t,h = 0,L) 
3 X W(t,h = 0,Ly 



(31) 



where means the n-th derivative with respect of h of the free energy density at h — (of course, x' 2 ' = — X, 
where \ is the susceptibility of the system). In the vicinity of the critical point this expression can be rewritten in 
the form 



B L (xx) 



d 4 X(xi,X2)/dx 



[d 2 X(x 1} x 2 )/dx 2 2 



(32) 



where X(x\,X2) is the finite-size scaling function of the free energy density. The exact form of this function follows 
from (fo) just by omitting the terms depending on j/oo in it, i.e. 



x{ Xl ,x 2 ) = --y iy£ 

1 xl 1 



3/2 



y L Li 2 (exp {-y/yl)) + Li 3 (exp (-y/yZ)) 



%vl 2 

From the above expression at the critical point it immediately follows that 

2 / q2 



B = B L {xi = 0) 



dyL 

0X2 



X1—X2—O 



d 2 y L 



X\— X2— 



(33) 



(34) 



By subsequent differentiation of the spherical filed equation for the finite system ( |18| ) it is easy to show that at the 
critical point dyL/dx 2 = 0, whereas 



d 2 y L 
dx\ 



16tt 



Xl =x 2 =o y 3 L /2 c coth ( A /yZ^/2) 



(35) 



Combining these results and having in mind that y^.c = 41n 2 r we obtain for the Binder's cumulant ratio at the 
critical point 



B 



2n 



\/5hi 3 



25.21657. 



(36) 



Having the exact solution for the spherical field equation and such a simple form for the free energy density, one can 
easily determine in an exact manner the behavior of other physically interesting quantities at T = T c . For example it 
is easy to show that the specific heat is of the form 
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cl(T c ) = ^-L- 1 ^Lk 2 c \ut, (37) 

and that the critical finite-size correlation length is (£l = Lj \fyZ [[52), ]37j]) 

UTc) = ^—L (38) 
2 mr 

(for explicit results of the behavior of ft under other geometries, boundary conditions and long-rangines of the 
spin-spin interactions see 02), ]37|, 138] , Q). 



V. CONCLUDING REMARKS 



In the present paper we present a hypothesis that in the vicinity of the bulk critical temperature T c of 0{n) 1 n > 2, 
systems with a film geometry L x oo d_1 the excess free energy (due to the finite size of the system) will be, under 
periodic boundary conditions, a monotonically increasing function of the temperature and the magnetic field if the 
finite system does not undergo a real phase transition of its own (i.e. when d < 3 for systems with short-range 
interactions) . As a function of the finite size L of the system the finite size scaling function of the excess free energy is 
expected to be monotonically increasing above and decreasing below T c . This hypothesis, together with the hypothesis 
that the Casimir force should be negative under periodic boundary conditions have been verified analytically on the 
example of the three-dimensional mean spherical model. It has been shown that the force is negative in the whole 
region of the thermodynamic parameters. In addition the universal Casimir amplitude A por and the Binder's cumulant 
ratio have been determined exactly in a simple close form and found to be A pcr = — 2£(3)/(57r) w —0.153051 and 
B = 2 7 r/( v / 51n 3 [(l + \/5)/2]) w 25.21657. For comparison we give the corresponding result for the Ising universality 
class, A pcr = -0.1526 ± 0.0010 ||, and B = 0.615 ± 0.003 @, (4l) obtained by Monte Carlo calculations. As 
we see, the value for the Casimir amplitude for the spherical model is surprisingly close (within the error bar) to 
the value reported above for the Ising model. The vast difference for the cumulant ratio indicates the lack of a real 
phase transition in the three dimensional spherical model film in comparison with the Ising like films. Actually, in 
three-dimensional Ising films the situation is more complicated |l2| . If the thickness of the film L is held constant and 
the other two linear dimensions D tend to infinity, the cumulant ratio converges to the two-dimensional Ising value 
(B = 0.615). However, if the ratio L/D is not too small, there exist crossover problems. In any case the value of B is 
between that one for the two-dimensional system and that one for the three-dimensional system (B — 0.47 |4j|). The 
value of B for the spherical model shows that the probability distribution at T c of the order parameter density is too 
different from a single Gaussian, where B — 0, or from a normalized sum of two Gaussians, where B = 2/3. This, of 
course, rises the question what then that distribution is, but this question is out of the scope of the current article. 
The situation reminds the one of Ising strips (no real phase transition in the system) with B — 2.46044 ± 0.00006 [[|, 
fio|| , The crossover problems in Binder's cumulant ratio can be studied within the spherical model, considering a 
3 + e dimensional film, e > (then in the finite system there will be a real phase transition) . This is also an interesting 
problem, especially if one takes into account that there are almost no exact results for the Binder's cumulant ratio, 
but it is again out of the scope of the current article. 

The results reported in the current investigation are in full agreement with the predictions of the finite-size scaling 
theory. Eqs.(|l7j), (|l8|), (19), ( |25| ) and ( |33| ) give the universal finite-size scaling function of the excess free energy, 
Casimir force and free energy density. It should be, however, emphasized that in contrast to the Ising-like case the 
excess free energy, and, therefore, the Casimir force in the absence of an external field tend to zero below T c not in an 
exponential in L way. For example, the finite-size scaling functions of the excess free energy and Casimir force tend 
to a constant below T c (see Eq. (31) in ||). The explanation of this behavior, which, we believe, is common for all 
0(n),n > 2 models, is based on the fact that due to the existence of soft modes in the system (spin waves) below T c 
and in the absence of an external field (h = 0) is identically infinite. If an external field is applied (h ^ 0) then 
£b < oo, and, of course, f cx — > again exponentially fast in L. 

Finally, it is worth mentioning the close parallel that exists between the properties of the C— function defined by 
Netto and Fradkin Q, see also Zabzin for a d-dimensional quantum system as a function of the temperature T 
and the properties of the excess free energy scaling function — A cx of the corresponding classical system as a function 
of L . If in the finite system a real phase transition does not exist, and if the system is somehow equivalent to the 
0(n),n > 2 system we have proposed some arguments that — X 0K is a monotonically increasing as a function of L^ 1 
above T c and decreasing below T c . We would expect the same to be true for the C-function of the corresponding 
quantum system as a function of T around its quantum critical point. If the classical system is equivalent to some 
Ising type model, the same type of arguments we have used for the 0(n),n > 2 models, taking into account the 
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lack of monotonicity of the correlation length in the vicinity of T c , lead to the hypothesis that — X cx will be a 
monotonic function of L^ 1 both below and above T c . For the corresponding C-function of a quantum system that 
has its mapping into a classical Ising system (according to the dimensional crossover rule) this means that C is a 
monotonically increasing function of the temperature both below and above its quantum critical point. This is indeed 
the case plotted in Fig. 2 in M for the quantum version of the two-dimensional Ising model. At the very end, we 
would like to stress that the relatively simple picture described here should probably change significantly, if the finite 
system undergoes a phase transition of its own. In that case the upper critical part of the excess free energy scaling 
function for 4 — e Ising model is known [po[ (up to a first order in e, e > 0). It shows a minimum in X CK , as a function 
of T slightly above T c . Unfortunately, no results are available for X cx when T < T c .l, where T c .l is the shifted critical 
temperature of the finite system. It is possible to investigate the above problems exactly within the spherical model 
with 3 + e infinite dimensions. We hope to return to this problem later. 
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Figure captions 

FIG.l. The universal finite-size scaling function of the excess free energy X°* as a function of the scaling variables 
x\ = L(K — K c )/K c and x-i = K c 1/>2 hL 5 ^ 2 . For a better visualization of the properties of X cx we have allowed h to 
change its sign. Of course, X cx is a symmetric function of x-i. 
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